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ABSTRACT 



I Context. The oxygen and nitrogen-to-oxygen abundances relation is characterized by a plateau around log(N/0) = -1.6 at metallicities 

lower than 12+log(0/H) ~ 8.25. However, there is a subset of Blue Compact Dwarf galaxies (BCDs) with unexpectedly high N/0 
values. 

^ ■ Aims. We want to study the spatial distribution of the physical properties and of oxygen and nitrogen abundances in three BCDs 

Q^l (HS 0128+2832, HS 0837+4717 and Mrk 930) with a reported excess of N/0 in order to investigate the nature of this excess and, 

particularly, if it is associted with Wolf-Rayet (WR) stars 

Methods. We have observed these BCDs by using PMAS integral field spectroscopy in the optical spectral range (3700 - 6900 A), 
mapping the brightest emission lines and, hence, their interstellar medium (ISM) physical-chemical properties (reddening, excitation, 
^ \ electron density and temperature, and O and N chemical abundances), using both the direct method and appropriate strong -line 

^ methods. We make a statistical analysis of the resulting distributions and we compare them with the integrated properties of the 

galaxies. 

Results. Outer parts of the three galaxies are placed on the "AGN-zone" of the [NIIJ/H vs. [OIII]/H diagnostic diagram most likely 
^ \ due to a high N/0 combined with the excitation structure in these regions. From the statistical analysis, it is assumed that a certain 

\^ , property can be considered as spatially homogeneous (or uniform) if a normal gaussian function fits its distribution in several regions 

of the galaxy. Moreover, a disagreement between the integrated properties and the mean values of the distribution usually appears 
1^ ' when a gaussian does not fit the corresponding distribution. We find that: 1) for HS 0128+2832, only the N/0 ratio as derived from 

\ the direct method is uniform; 2) for HS 0837+4717, all properties are uniform; and 3) for Mrk 930, the uniformity is found for all 

• , parameters, except for electron density and reddening. The rotation curve together with the Ha map and UV images, reveal a perturbed 

morphology and possible interacting processes in Mrk930. Finally, if the N/0 is constant at spatial scales of the order of several kpc, 
' as our analysis indicates, we conclude that the number of WR stars estimated from spectroscopy, is not sufficient to pollute the ISM 

and to produce the observed N/0 excess in these objects. 

Key words, galaxies : evolution - galaxies : abundances - galaxies : starbursts - galaxies : kinematics and dynamics - galaxies : 
stellar content 

X 

■ 1, Introduction characterized very accurately by analyzing their optical spectra 

, r , r ' 11. r 1 in a wide sample of these objects (e.g. Terlevich et al., 1981: 

The understanding of the formation and evolution of staxburst ^^^^^^ ^ q^^^^ 2000). Since the metalhcity distribution of these 
galaxies has been based mainly so far on the study of differ- ^^-^^^^ ^ ^^^^^ ^^^^^^ . ^^^^ ^^^j^ ^^^^ 

ent integrated properties of their star-forming knots, which de- -^^^^^ ^^^^^ 1^^^ ^^^^-^^^ ^^^^ metal-rich giant Hn 
pend on their content of stars gas, and dust. In the optical spec- -^^^ galaxy-disks and consequently their electron temper- 
trum of these knots, the blue light emitted by the massive young ^^^^^^ ^^^^ ^- ^^^ ^^-^ ^^^^ ^j^^ determination of chemi- 
star clusters and the emission lines coming from the ionized gas abundances using the so-called direct method (i.e. using col- 
stand out. These lines provide information about the physical Usionally excited lines and electron temperatures) much more 
properties and the chemical abundances of the gas and the stars. ^^^^^^^^ ^^^^ ^^^^^ metallicity regimes (e.g. Perez-Montero 
Between the different types of starburst galaxies. Blue Compact ^ 2003, Hagele et al., 2006). This feature of BCDs, to- 
Dwarf galaxies (BCDs), a lso known as Hii galaxies, have been g^^j^^^ ^-^^^ ^^^^-^ ^^^^^^^ ^^p^^^ ^-^^^^^^^^ ^^^^ ^j^^^ 

Send offprint requests to: E. Perez-Montero ideal scenarios to study the star formation processes in the local 

* CONICET Argentina universe and m low metallicity environments resembling those 
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in a younger Universe. Particularly, it is interesting to look into 
the ionic abundance ratios between species with an assumed dif- 
ferent stellar origin to probe the chemical evolution of these ob- 
jects. 

This is the case of the nitrogen-to-oxygen ratio (N/O), since 
their stellar yields are quite different and mostly proceeding from 
low- and intermediate-mass and massive stars, respectively, and, 
therefore, it gives important clues about the star formation rate 
and history of star-forming galaxies (Molla et al., 2006). This 
ratio has been well studied since some decades ago. Edmunds & 
Pagel (1978) and AUoin et al. (1979) concluded that the constant 
value of N/O [log(N/0) ~ -1.6] in the low metallicity regime (i.e. 
for 12-Flog(0/H) < 8) is consistent with a primary origin of the 
nitrogen, while the correlation between N/O and 0/H at higher 
oxygen abundances corresponds to a secondary production of 
N. These findings were later confirmed by McCall, Rybsky & 
Shields (1985) and Vila-Costas & Edmunds (1993), among oth- 
ers. 

Traditionally, both the plateau and the dispersion found in 
the N/O vs. 0/H plot at the low metallicity regime have been ex- 
plained in terms of the delay between the ejections of O and N, 
which appears if O is produced by massive stars and N by low- 
and intermediate-mass stars. In that case, different star forma- 
tion histories included in chemical evolution models may change 
the exact value of 0/H for which the secondary N would appear 
and, therefore, the observed dispersion is easily explained, such 
as Henry et al. (2000), Prantzos (2003) and Molla et al (2006) 
found. At the same time the plateau appears as a consequence of 
a very low and continous star formation rate and of metallicity- 
dependent stellar yields which, as expected, give a high propor- 
tion of primary N when 0/H is low, and a higher secondary N 
production when 0/H increases. 

Besides, the relative production of N is sensitive to the as- 
sumed theoretical stellar yields. Thus, Renzini & Voli (1981) 
predicted a certain amount of primary N produced by the low- 
and intermediate-mass stars. Other works also give a propor- 
tion of primary N for this stellar mass range, such as Van der 
Hoek & Groenewegen (1997) and, more recently, Gavilan et al. 
(2006). These stellar yields, computed for a wide range of metal- 
licities and masses, and by simplified synthesis stellar models 
have been also corroborated by precise calculations from the 
stellar evolution field (Marigo, 2001; Ventura et al., 2002;, Dray, 
2003). Additionally, Meynet & Maeder (2002) and Chiappini et 
al. (2006) predict an enhancement of the N production in the low 
metallicity regime assuming rapid rotation of massive and inter- 
mediate mass stars, fitting the observations coming from the halo 
stars of our Galaxy. 

Moreover, in any case, by assuming a recursive star forma- 
tion history with several bursts is also possible to find some 
extra-enhancements of the N/O ratio in the starburst phases: 
O increases when massive stars die, decreasing simultaneously 
N/O, and then, when intermediate mass stars die, N increases 
again. Hence, some peaks should appear in the 0/H vs. N/O di- 
agram (Gamett, 1990; Pilyugin, 1992). Several authors have ap- 
pealed to hydro-dynamical effects, such as outflows of enriched 
gas by super-massive galactic winds (van Zee et al., 1998) or in- 
flows of metal-poor gas (Koppen & Hensler, 2005), which makes 
the metallicity of the galaxy to be lower, keeping the N/O ratio 
at a high value, to explain the observations. 

In this context, there is a subsample of BCDs which present 
high values of N in comparison to the expected value for their 
metal content, even taking into account the observational disper- 
sion (Perez-Montero & Contini, 2009). Two examples are Mrk 
996 [Thuan et al., 1996, 12-Flog(0/H)=7.96, log(N/O)=-0.80] 



and NGC 5253 [Walsh & Roy, 1987, 1989; Kobulnicky et al., 
1997, 12-Flog(0/H)=8.12, log(N/O)=-0.83 in the region HII-1]. 
More recently, "Green pea" galaxies, which resemble in many 
aspects BCDs, have been found to be also depicted by low metal- 
licities and high N/O ratios (Amorm, Perez-Montero & Vilchez, 
2010). In a number of these metal-poor BCDs with an enhanced 
N/O ratio has been detected broad prominent emissions at ~ 
4650 A and ~ 5808 A. These broad features are emitted by Wolf- 
Ray et (WR) stars, being characterized by the ejection of material 
via strong stellar winds which pollute the interstellar medium 
(ISM) with the products of the H and He-burning in the outer 
layers of these stars. This suggests that WR stars could be also 
the responsible for the high N abundances in these galaxies. 

Among BCDs showing high N/O and where WR stars were 
detected, are UM420, UM448, the merger compact group Mrk 
1089 (all these in Guseva et al., 2000), HS 0837-h4717 (Pustilnik 
et al., 2004) and NGC 5253 (Lopez-Sanchez et al., 2007; 
Monreal-Ibero et al., 2010). In the study of galaxies with WR 
signatures from the Sloan Digital Sky Survey (SDSS) carried out 
by Brinchmann et al. (2008), it is also found that WR galaxies 
show an elevated N/O relative to non-WR galaxies. Additionally, 
a N/O excess is reported in other six BCDs with log(N/0) > - 
1.3 and a prominent WR blue bump in their optical integrated 
spectra (Hagele et al., 2006, 2008, Perez-Montero et al., 2010). 
On the other hand, the high N/O ratio found in metal-poor halo 
stars indicates that another mechanism polluting the ISM with 
N, other than WR star winds, might exist. 

Nevertheless, most of the contributions which establish the 
relation between the local N pollution and the presence of WR 
stars are based on long-slit integrated observations and, there- 
fore, almost no information about the relative spatial position 
of massive stars and possible chemical inhomogeneities exist in 
this sample of BCDs. One of the few examples is described in 
Kehrig et al. (2008) who nd an excess of N/O simultaneously 
with the detection of WR stellar population across the BCD 
IIZw70 by means of integral eld spectroscopy (IFS). These are 
also the cases of NGC5253 (Monreal-Ibero et al., 2010) and 
ICIO (Lopez-Sanchez et al., 2010). IFS constitutes a powerful 
tool to study in detail the spatial distribution of the different 
physical properties and chemical abundances across the objects. 
Besides, an appropriate statistical analysis of these distributions, 
helps to distinguish between the spatial uniformity of some prop- 
erties and to detect local inhomogeneity or pollution of chemical 
elements in the ISM of galaxies, and possibly to link them with 
the positions and properties of the star-forming knots. 

In this work, we use IFS to analyze three BCDs (HS 
0128-h2832, HS 0837-h4717, and Mrk 930) selected because 
of their overabundance of N reported in the literature. In HS 
0837-^4717 and Mrk 930, the presence of WR stars has been 
confirmed by long-slit spectroscopy. IFS lead to a 2D spa- 
tial characterization of both the stellar population and the ISM 
physical-chemical properties allowing us to check whether the 
N/O excess is spatially correlated with the WR stars or if it is 
due to a more global process affecting the whole galaxy. 

The paper is organized as follows. In Section 2, we describe 
our sample of objects as well as the observations and data re- 
duction. In Sections 3 and 4, we present and discuss our re- 
sults. Finally, Section 5 summarizes the main conclusions de- 
rived from this work. 
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Table 1. List of the observed objects, with some additional basic information taken from Izotov & Thuan. (2004; HS 0128+2832), 
Pustilnik et al. (2004; HS 0837+4717) and Izotov & Thuan (1998; Mrk 930). The table includes for each object coordinate position, 
adopted distance, redshift, logarithm of Ha flux in erg/s/cm^A, and oxygen and nitrogen chemical abundances. These abundances 
have been recalculated using the procedure described in the text and taking the emission-line data provided in the same sources. 



Object ID 


Other names 


R.A. (2000) 


6 (2000) 


D(Mpc) 


redshift 


log F(HQf) 


12+log(0/H) 


log(N/0) 


HS 0128+2832 




01h31m21.3s 


+28d48ml2s 


63.1 


0.016130 


-13.23 


8.12 + 0.05 


-1.23 + 0.14 


HS 0837+4717 


SHOC 220 


08h40m29.9s 


+47d07ml0s 


166.0 


0.041950 


-13.06 


7.62 + 0.06 


-0.75 + 0.20 


Mrk 930 




23h31m58.3s 


+28d56m50s 


83.2 


0.018296 


-12.71 


8.07 + 0.03 


-1.42 + 0.10 



2. Data 

2.1. The sample 

Three BCDs with difl'erent degrees of N overabundance as de- 
rived using integrated long-slit spectroscopy available in the lit- 
erature and visible from the Northern hemisphere at the epoch 
of observations were selected. The main properties of these ob- 
jects as taken from the literature are listed in Table [T] includ- 
ing names, position, adopted distance, redshift and Ha flux in 
the slit without aperture correction. The available data comprise 
the required emission lines to derive at least one electron tem- 
perature allowing the direct and precise determination of the O 
and N chemical abundances. These were recalculated using the 
methodology described in the section 3.6, and they are also listed 
in Table [TJ In Fig.[Tl we show the derived O abundance and N/0 
ratio for these three objects compared to the values from the star- 
forming galaxies of the SDSS (Amorm et al., 2010). As we see, 
all of them present values of N/0 higher than the average N/0 
for their same metallicity regime. 

HS 0128+2832 is a compact galaxy studied by Izotov & 
Thuan (2004), who derived a low O abundance [12+log(0/H) 
= 8.12] combined with an elevated abundance of N [log(N/0) = 
-1.23]. No trace of WR stars was reported in the integrated spec- 
trum of this object. HS 0837+4717 is apparently another com- 
pact galaxy which has been extensively studied by Pustilnik et al. 
(2004), who found very low O abundance [12+log(0/H) = 7.65] 
and a very high value of the relative abundance of N [log(N/0) 
= -0.75]. These authors detected the WR blue bump which was 
used to estimate approximately 1000 WR stars in the starburst. 
Although the stellar origin of the Hen line shown by Pustilnik 
is not clear, this galaxy also appears in the catalog of WR ob- 
jects in the SDSS (WRIOI) as a class 2 object (i.e. convincing 
WR feature but not obvious after continuum subtraction). Our 
third object, Mrk 930, presents a quite disturbed morphology, 
characterized by an ellongated shape and whose main burst of 
star formation is in the south, but other several fainter knots are 
also well visible in the north part. Fig. [2] shows an ACS-HSTQ 
image taken with the F140LP filter centered in the ultraviolet at 
a wavelength ~ 1400 A, where it is remarkable the presence of 
these super star clusters along the galaxy in a complex structure. 
Guseva et al. (2000) analyzed both the blue and red WR bumps 
detected by Izotov & Thuan (1998) for this galaxy and they es- 
timated around 2000 WR stars, although this is very uncertain 
due to the noise of the spectrum. A very noisy hint of this WR 
feature is also pointed out by Adamo et al. (201 1) using long-slit 
spectroscopy of the brightest knot of this galaxy. The chemical 
analysis carried out by Izotov & Thuan (1998) for this galaxy 
revealed a low metallicity [12+log(0/H) = 8.07] and a value of 
the N/0 ratio slightly enhanced [log(N/0) = -1.42]. 



^ This image has been taken from the Multimission Archive at STScI 
(MAST) in the webpage http://archive.stsci. edu/, 




12+log(0/H) 



Fig. 1. Relation between O abundance and N/0 ratio as derived 
from emission-line spectra of star-forming objects of the SDSS. 
The red squares represent the objects of the sample described in 
the text, with their values recalculated using the emission lines 
reported in the literature 

The study of these three galaxies by means of IFS constitutes 
an ideal test to find out to what extent the relative enrichment of 
N found in some BCDs is related to the presence of WR stars in 
their knots of star formation. 

2.2. Observations 

Observations were carried out using the integral field unit 
(IFU) Postdam Multi- Aperture Spectrophotometer (PMAS), de- 
veloped at the Astrophysikalisches Institut Potsdam (Roth et al. 
2005). PMAS is attached to the 3.5 m Telescope in the Calar 
Alto Observatory (Almeria, Spain). The data were acquired on 
2008, October, 29th, under good conditions of transparency and 
with a typical seeing close to T\ To avoid major diff'erential at- 
mospheric reff'raction eff'ects, all expositions were taken at an air 
mass lower than 1.2. 

One single pointing was taken for each galaxy, covering in 
all cases the most intense burst of star formation and its sur- 
roundings. 

The log of observations is given in Table [21 The V600 grat- 
ing, with a dispersion of 1.6 A/pixel in the 2x2 binning mode 
to reduce the reading time and readout noise, was used in two 
spectral ranges, covering both the blue (3700-5200 A) and red 
(5350-6850 A) sides. The PMAS spectrograph is equipped with 
256 fibers coupled to a 16 x 16 lens array. Each fiber has a spa- 
tial sampling of f x on the sky, resulting in a field of view 
of 16'' X 16''. This field is shown in relation to the size of Mrk 
930 in Fig.[2j At the reported distances of the three objects (see 
Table 1), each spaxel size of 1" corresponds to 330 pc in HS 
0128+2832, 840 pc in HS 0837+4717, and 380 pc in Mrk 930. 
Calibration images were taken following the science exposures 
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RA 02000) 



Fig. 2. HST-ACS image of Mrk 930 taken with the F140-LP fil- 
ter. The square shows the field of view of the PMAS instrument 
encompassing the whole galaxy. The figure also shows the com- 
plex structure of the ionizing stellar clusters in this galaxy. North 
points towards up and east to the left. 

Table 2. Log of observations. 



Object 


B exptime 


R exptime 




(s) 


(s) 


HS 0128+2832 


2x1200 + 3x900 


2x1200 + 3x900 


HS 0837+4717 


2x1200 + 1x900 


2x1200 + 1x900 


Mrk 930 


2x1200 


2x1200 



and consisted of emission line lamp spectra (HgNe), and spectra 
of a continuum lamp needed to locate the 256 individual spectra 
on the CCD. Observations of the spectrophotometric standard 
stars BD +28°4211 and Hz44 were taken during the observing 
night for flux calibration. 



2.3. Data reduction 

We reduced the data using the software R3D (Sanchez 2006). 
Difl'erent exposures taken at the same pointing were combined 
using IRAI0 tasks. The expected locations of the spectra were 
traced on a continuum-lamp images taken before each target ex- 
posure. After bias subtraction, we extracted the target spectra 
by adding the signal from the 5 pixels around the central traced 
pixel (which is the total object spectrum width). With exposures 
of Hg and Ne lamps taken immediately after the science expo- 
sures, the spectra were wavelength calibrated. We checked the 
accuracy of the wavelength calibration using sky emission lines, 
and find typical deviations of 1 A and 0.5 A for the blue and red 
ranges, respectively. The eff'ective spectral resolution, derived by 
measuring the width of the arc lines is 3.5 A FWHM (corre- 
sponding to an instrumental dispersion about 90 km/s at 5007 A 
and 70 km/s at Ha wavelength). The continuum-lamp exposure 
was also used to determine the response of the instrument for 



^ IRAF is distributed by the National Optical Astronomy 
Observatories. 



each fiber and wavelength (the so-called flat- spectra). We used 
these flat- spectra in order to homogenize the response of all the 
fibers. For the standard star observations we co-added the spectra 
of the central fibers and compared the one-dimensional standard 
star spectrum with table values to create a sensitivity function. 
The spectra were flux calibrated using IRAF. 

The reduced spectra were contained in a data cube for each 
object and spectral range and they were sky- subtracted and cor- 
rected for the efl'ect of DAR using the R3D package (Sanchez 
2006). 

3. Results 

3.1. Line measurement and Ha maps 

We fitted each emission-line profile on the extracted one- 
dimensional spectra of each spaxel using a gaussian function 
over the local position of the continuum. 

We used an automatic routine based on the IRAF task splot 
to measure the flux of most of the emission-lines. For the low 
S/N emission-lines (e.g. [Oiii] 4363 A) and those lying close to 
the border of the CCD (e.g. [On] 3727 A for HS 0128+2832 
and Mrk 930), the fitting procedure was repeated using a manual 
measurement. 

The measured flux of the B aimer emission lines is some- 
times underestimated because of the the presence of an under- 
lying stellar continuum (e.g. Diaz, 1988). For objects showing 
high EW(HyS) values (223 A in HS 0128+2832, 230 A in HS 
0837+4717, and 93 A in Mrk 930), the eff'ect of the steflar ab- 
sorption is negligible due to the low typical values of 1-2 A for 
the hydrogen line EWs in absoprtion (see e.g. McCall et al. 
1985). However, to be careful, we carried out an eye-inspection, 
fiber-by-fiber, of the HJS and Hy emission lines, and we did not 
find any apparent stellar absorption feature (e.g. wings of ab- 
sorption lines) underlying these emission lines. 

We calculated the statistical error of the line fluxes, cr/, using 
the expression cr/ = (TcN^/^[1+EW/NA]^/^ (see Perez-Montero 
& Diaz, 2003) where cTc represents a standard deviation in a box 
centered close to the measured emission line, N is the number 
of pixels used in the measurement of the line flux, EW is the 
equivalent width of the line, and A is the wavelength dispersion 
in A/pixel. This expression takes into account the error in the 
continuum and the photon counts statistics of the emission line. 
The error measurements were performed on the extracted one- 
dimensional spectra. These associated errors are quite difl'erent 
depending on the observed object. Hence, the average relative 
error of HyS is 10% in HS 0128+2832 and HS 0837+4717, but 
only of 5% in Mrk 930. Besides, to minimize errors in the ratios 
between a certain emission line and UJ3, we always take first its 
ratio in relation to the closest hydrogen emission line (i.e. Ha in 
the case of [Nii] and [Sii]) and then we renormalize it using the 
corresponding theoretical B aimer ratio (i.e. at the electron tem- 
perature reported in the long-slit observations for each object). 
We have checked that the variation of this temperature across the 
field of view of the instrument does not introduce errors in the 
theoretical ratio higher than those associated with the flux of the 
emission lines. 

Ha emission line maps (continuum subtracted and extinction 
corrected) are shown in Fig. [3] In the same figure, the contours 
in grey solid line show the isophotes in units of 12.5 % times 
the emission of the maximum of Ha in each galaxy. In the panel 
of Mrk 930, we also show the contours extracted from the HST- 
ACS image shown in Fig. 2. Although the combined low spatial 
resolution of the instrument (1 "/spaxel) and poor seeing during 



4 



E. Perez-Montero et al: A PMAS study of high N/0 BCDs 




HS 0837+471 7 



2E-15 
lE-15 





Fig. 3. Extinction corrected Ua maps of HS 0128+2832, HS 0837+47 17, and Mrk 930 from left to right respectively. The grey solid 
line contours show the isophotes in units of 12.5 % of the maximum of Ha. In the Mrk 930 image, black solid lines represent the 
contours extracted from the ACS-HST UV image shown in Eig. |2l which represent the position of the brightest super star clusters. 
In all images, each spaxel has 1" of resolution, north points to up and east to left. Units are in erg / s / cm^ / A. 



HS 0837+4717 



i 



0.2 0.4 0.6 0.8 

cm 




i; 





Fig. 4. Reddening constant maps of HS 0128+2832, HS 0837+4717, and Mrk 930 from left to right, respectively. The solid line 
contours are those described in Eig. 3. The histogram plots above represent the distribution of the extinction in different regions: 
Dark blue bars represent the spaxels in Region 1 (spaxels with a E(HQf) larger than a 25 % of the maximum), light blue bars, in 
Region 2 (spaxels with E(Ha) larger than 12.5 % of the maximum), and white bars. Region 3 (all the spaxels in the observed field 
with enough S/N). 



the observing night (~ 2") do not allow a precise analysis of the 
morphology of the galaxies, it is seen that in HS 0128+2832 
and HS 0837+4717, the central compact burst dominates the 
emission of the galaxy, while Mrk 930 displays a more extended 
distribution of the star formation in agreement with what is ob- 
served in the UV image (see Eig. [2]). 



3.2. Statistical distributions and integrated spectra 

Line intensities with their corresponding errors measured for 
each spaxel have been used to obtain the maps of several phys- 
ical properties and chemical abundances. Only those emission 
lines with a signal-to-noise ratio (S/N) > 5 have been consid- 
ered. In order to extract statistical information about the be- 
haviour of the various derived physical properties and chemi- 
cal abundances, we also provide their distribution histograms, 
taking into account diff'erent regions around the Ha peak in- 
tensity, F(}la)peak- Hence, hereafter in all histograms, dark blue 
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Table 3. Results of the Lilliefors test and the gaussian fits to the histograms of the derived physical conditions and chemical 
abundances in the three studied objects for diff'erent subsets of spaxels: Region 1: spaxels with a relative F(HQf) larger than 25 % of 
the maximum ¥(tla) in the galaxy, Region 2: larger than 12.5% of the maximum flux, Region 3: all spaxels in the field with enough 
signal-to-noise for the corresponding involved emission lines. For each region the number of spaxels considered for the statistics 
is done. Sign, stands for the significance level of the null hypothesis in the Lilliefors test. If the null hypothesis cannot be rejected 
with more than a 10 % of confidence, then the gaussian fit is assumed and the corresponding mean and dispersion values are taken 
(in bold font). Otherwise, the mean and dispersion are those of the not fitted distribution. 



Object 




HS 0128+2832 


HS 0837+4717 




Mrk 930 




Region 




1 


2 


3 


1 


2 


3 


1 


2 


3 


log F(HQf) (erg/s/cm^/A) 


-i J.UU 


-iZ.OJ 


1 9 ^7 
-iZ.J / 


1 9 8zL 


1 9 


1 9 ZL9 


- 12.06 


-11.95 


-11.87 


C(Hyg) 


Spaxels 


1 

1 J 






1 8 


H-VJ 


77 
/ / 


42 


73 


157 




Mean 


yj.D 1 


40 






7Q 


S7 


0.48 


0.46 


0.51 




Dispersion 


U.ZJ 


U.ZO 


n 9Q 








0.24 


0.24 


0.31 




Sign. (%) 


< iU 


< iU 


< iU 


71 


oU 


40 


< 10 


< 10 


< 10 


log ([Oii]/[Oiii]) 


Spaxels 


14 


9Q 




1 7 


98 


D 1 


47 


66 


143 




Mean 


n Q7 




-U.oU 


n QzL 


-U.oo 


-U.o,3 


-0.31 


-0.33 


-0.22 




Dispersion 


U. iU 


n 9n 
u.zu 




91 
U.Zi 




U.ZO 


0.18 


0.18 


0.21 




Sign. (%) 


< iU 


< iU 


< iU 


< iU 


1 1 

1,3 




38 


64 


< 10 


n([Sii]) (cm-^) 


Spaxels 














38 


67 


107 




Mean 














101 


213 


271 




Dispersion 














135 


275 


296 




Sign. (%) 














< 10 


< 10 


< 10 


t([Oiii]) (lO^K) 


Spaxels 


1 


99 
ZZ 


99 
ZZ 


iU 


1 9 
iZ 


1 


39 


52 


52 




Mean 


1 '^^ 
i.Jj 


i. J / 


1 '^7 
i. J / 


1 /i/i 


1 /id 




1.26 


1.29 


1.29 




Dispersion 


U.iU 


n 1 Q 


n 1 Q 








0.18 


0.21 


0.21 




Sign. (%) 


< iU 


^ 1 n 
< iU 


< iU 


LL 


51 




11 


15 


15 


12+log(0/H) 


Spaxels 


Q 

y 


91 


91 


1 


1 9 


1 


35 


47 


47 


(direct method) 


Mean 


8.00 


7.99 


7.99 


7.83 


7.82 


7.85 


8.18 


8.15 


8.15 




Dispersion 


0.09 


0.10 


0.10 


0.15 


0.14 


0.18 


0.16 


0.18 


0.18 




Sign. (%) 


< 10 


< 10 


< 10 


41 


15 


11 


55 


32 


32 


12+log(0/H) 


Spaxels 


14 


29 


50 


16 


27 


32 


41 


79 


79 


(R23-P method) 


Mean 


8.01 


7.99 


7.96 


7.88 


7.90 


7.92 


8.18 


8.14 


8.14 




Dispersion 


0.09 


0.11 


0.28 


0.18 


0.19 


0.23 


0.16 


0.19 


0.19 




Sign. (%) 


< 10 


25 


< 10 


98 


17 


18 


13 


< 10 


< 10 


log(N/0) 


Spaxels 


13 


17 


17 


9 


11 


12 


36 


49 


51 


(direct method) 


Mean 


-1.19 


-1.17 


-1.17 


-0.82 


-0.81 


-0.82 


-1.43 


-1.44 


-1.43 




Dispersion 


0.11 


0.16 


0.16 


0.16 


0.15 


0.14 


0.13 


0.19 


0.18 




Sign. (%) 


72 


67 


67 


82 


83 


57 


22 


12 


38 


log(N/0) 


Spaxels 


15 


25 


25 


17 


28 


34 


47 


84 


124 


(N202 method) 


Mean 


-1.18 


-1.12 


-1.12 


-0.85 


-0.80 


-0.79 


-1.36 


-1.43 


-1.41 




Dispersion 


0.25 


0.27 


0.27 


0.26 


0.28 


0.29 


0.24 


0.17 


0.11 




Sign. (%) 


< 10 


< 10 


< 10 


< 10 


18 


27 


28 


< 10 


< 10 



bars represent the spaxels for which F(HQf) > 25% • ¥(}\a)peak 
and light blue bars, those where F(HQf) > 12.5% • ¥(]:{a)peak- 
Finally, white bars represent the distribution of the correspond- 
ing physical-chemical property in all spaxels with enough S/N 
over the all IFU area. In Table [3] we show the total Ha 
extinction-corrected fluxes in logarithm units, as measured in 
each of these regions for the three studied galaxies. The ratios 
between the Ha fluxes measured in the whole field of view of the 
IFU and the fluxes measured in long-slit observations, listed in 
Table 1, can be considered as aperture correction factors, giving 
an estimate for the Ha flux lost when using long-slit observa- 
tions. In our cases, this gives ratios of 4.6 in HS 0128+2832, 4.4 
in HS 0837+4717, and 6.9 in Mrk 930. In the case of Mrk 930, 
this factor is almost three times the aperture correction factor 
estimated by Guseva et al. (2000). 

In order to ease the comparison between the studied distribu- 
tions and the position of each defined region, we have overplot- 
ted in all the maps the corresponding Ha galaxy contours, with 
the exception of Mrk 930, for which we have used the contours 
from the ACS image, with better spatial resolution, easing the 
identification of the individual knots. 



In order to ascertain if a physical-chemical property can be 
considered as homogeneous across the IFU area, we have fol- 
lowed a statistical criterion. We have plotted the spaxel distribu- 
tion of a given property in each region of each galaxy and we 
have studied if that distribution can be fitted by a normal func- 
tion. As a first-order approximation, it can be assumed that the 
measured variations around the mean value of the distribution 
have mainly a statistical origin and, hence, a normal distribution 
of measurements/derivations can be considered. With this aim, a 
Lilliefors test (Lilliefors, 1967) was carried out for each of the 
distributions of the derived properties and abundances in each 
of the regions defined above. For this test, the null hypothesis 
is that the data come from an unspecified normal distribution. 
We then assumed that if the null hypothesis cannot be rejected 
at the 10% of significance level, the distribution is considered as 
normal. In that case, the fitted gaussian is plotted in the corre- 
sponding histogram. The number of spaxels used to do this test 
in each region of the IFU, with the mean value, the dispersion of 
the fitted gaussian and the significance level are listed in Table [3] 
and the corresponding values are noted in bold font. Otherwise 
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Table 4. Dereddened emission line intensities, relative to I(HyS) = 1000, and physical properties and oxygen and 
nitrogen chemical abundances with their corresponding propagated errors, as derived from the integrated spectra 
of the studied galaxies. Mrk 930 has been divided between the brightest southern knot, and those situated in the 
northern region. 





HS 0128+2832 


HS 0837+4717 


Mrk 930-S 


Mrk 930-N 


^(A) 


I(^) 


I(^) 


I(^) 


1(A) 


3727 [On] 


788 + 116 


1045 + 87 


2672 + 43 


2532 + 58 


4102 U6 


270 + 9 


271 + 22 


277 + 10 


292 + 18 


4340 Hy 


455 + 11 


470 + 24 


484 + 11 


474 + 18 


4363 [Oiii] 


100 + 8 


149 + 12 


54 + 6 


40 + 5 


4861 H/? 


1000 + 27 


1000 + 45 


1000 + 15 


1000 + 20 


4959 [Oiii] 


2269 + 48 


2021 + 77 


1593 + 22 


1452 + 21 


5007 [Oiii] 


6962 + 139 


6064 + 195 


4817 + 59 


4441 + 65 


6563 Ha 


2811 +66 


2774 + 102 


2830 + 47 


2849 + 53 


6584 [Nil] 


78 + 8 


152+ 10 


128 + 10 


136 + 20 


6717 [Sii] 


231 + 20 


- 


257 + 10 


309 + 30 


6731 [Sii] 


176 + 20 


- 


198 + 10 


215 + 20 




0.14 + 0.05 


0.80 + 0.10 


0.43 + 0.07 


0.32 + 0.08 


log([Oii]/[Oiii] 


-1.07 + 0.09 


-0.89 + 0.07 


-0.37 + 0.02 


-0.38 + 0.01 


n([Sii]) fcm~^ ) 


200+?SX 




< 190 


32+?^ 

^^-19 


t([Oiii]) (10^ K) 


1.32 + 0.05 


1.69 + 0.10 


1.20 + 0.05 


1.11 +0.05 


t([Oii]) (10^ K]^ 


1.12 + 0.08 


1.45 + 0.20 


1.16 + 0.05 


1.11 +0.07 


t([Nii]) (10^ Kf 


1.25 + 0.03 


1.41 +0.04 


1.19 + 0.03 


1.14 + 0.03 


12+log(0+/H+) 


7.09 + 0.24 


6.90 + 0.21 


7.69 + 0.10 


7.89 + 0.13 


12+log(02+/H+) 


7.99 + 0.06 


7.67 + 0.07 


7.95 + 0.06 


8.01 + 0.07 


12+log(0/H) 


8.04 + 0.08 


7.73 + 0.10 


8.14 + 0.05 


8.24 + 0.12 


12+log(0/H) (R23-P) 


7.95 


7.90 


8.19 


8.10 


12+log(N+/H+) 


5.95 + 0.10 


6.14 + 0.05 


6.14 + 0.03 


6.22 + 0.03 


log(N/0) 


-1.14 + 0.25 


-0.76 + 0.23 


-1.47 + 0.10 


-1.59 + 0.15 


log(N/0) (N202) 


-1.14 


-0.98 


-1.41 


-1.50 



^ Derived from t([Oiii]) using the models from Perez-Montero & Diaz (2003) and Perez-Montero & Contini (2009) for 
t([Oii]) and t([Nii]), respectively. 



(i.e. if the significance level is lower than 10%), the mean and 
dispersion values are those of the not fitted distribution. 

A comparative analysis between the distribution of the de- 
rived properties in the spaxels and their values in the integrated 
spectra has been done by means of the co-addition of the emis- 
sion of those spaxels located in the central region (i.e. presenting 
Ha emission higher than a 25 % of the intensity in the peak of 
each galaxy). Contrary to the statistical distribution method de- 
scribed above, this procedure gives a flux-weighted average of 
the properties of the galaxy and it is more consistent to estab- 
lish a comparison with long-slit spectra. In the case of ]V[rk930, 
we have extracted two spectra: one around the southern brightest 
knot and a second one close to the northern knots. The relative 
dereddened emission lines and the derived properties for each of 
the integrated spectra, as described in the subsections below are 
summarised in Table IH 



3.3. Reddening correction and c(Hf3) maps 

For each fiber spectrum we derived its corresponding reddening 
coefficient, c(H/5), using the value of the B aimer decrement de- 
rived from Ha^/HyS, Hy/HyS, and H^/HyS flux emission-line ratios, 
as compared to the theoretical values expected for recombination 
case B from Storey & Hummer (1995) at the electron density 
and temperature reported in the literature using long-slit obser- 
vations, and applying the extinction law given by Cardelli et al. 
(1989) with Ry = 3.1. Thus, the fluxes of the emission lines for 
each fiber were corrected for extinction using their correspond- 
ing c(Uj3) values. 



The distribution functions and the maps of the reddening 
constants, c(HyS), are shown in Fig. [4] along with the same Ha 
and UV (in the case of IMrk 930) contours described above. The 
typical errors of the derived reddening constants are 0.2 in the 
three objects. A gaussian fits the distribution of the three consid- 
ered regions of HS 0837+4717 with a mean value quite similar 
to that found in the integrated spectrum [c(Hy5) = 0.80]. On the 
contrary, in HS 0128+2832 none of the three distributions are 
fitted by a gaussian function and the derived value for the inte- 
grated spectrum of the inner-most region is much lower (0.14) 
than the mean of the distribution for the same region (0.37), 
what could be indicative of a lower extinction in the Ha peak 
region. This behaviour also appears in the histogram of Region 
1 for this galaxy, where a double-peak structure can be seen. In 
IMrk 930 no gaussian fits any of the three distributions, which 
are centeres at low values of the reddening, although can be seen 
some spaxels with reddening c(Uj3) > 1.0 breaking the normal 
distribution in all the regions. By visual inspection of the cor- 
responding c(HyS) map, these appear above all in the south-west 
part of the galaxy, resembling somehow the ring-like structure 
shown by Adamo et al. (201 1). 

3.4. Ionization structure and excitation 

Diagnostic diagrams (BPT, Baldwin, Philips & Terlevich, 1981) 
are frequently used to distinguish between diff'erent ionization 
sources of the surrounding gas emitting the detected bright emis- 
sion lines. Active galactic nuclei (AGNs), shocks or massive 
stars can be selected, by comparing the ratios of low excitation- 
to-Balmer emission lines, like [Nii] at 6584 A, [Sii] at 6717 A, 
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Fig. 5. Relation between log([Nii]/HQf) and log([Oiii]/H/5). one of the so-called BPT diagrams for the individual spaxels of the 
three studied galaxies. Dark blue squares represent the spaxels in Region 1, the light blue, Region 2, and the white all the rest. 
Large stripped squares represent the result of the integrated spectrum of Region 1. The red dashed line represents the empirical 
curve defined by Kauff'man et al. (2003) to separate star-forming galaxies from active galactic nuclei. The solid line represents the 
theoretical separation defined by Kewley et al. (2001). 
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Fig. 6. Distribution histograms (top) and maps (bottom) of the excitation maps in terms of log([Oii]/[Oiii]) in HS 0128-^2832, HS 
0837-h4717, and Mrk 930 from left to right respectively. The solid lines show the same contours as in Fig. 3. Units are in dex. The 
distribution regions are the same defined in Fig. |4] 



6731 A or [Oi] at 6300 A, to Ha in relation to a high excitation- 
to-Balmer emission line ratios, like [Oiii] at 5007 A to H/S. In 
all studies about our three galaxy sample, their integrated spec- 
tra show a clear star-forming ionization origin in all these three 
diagrams. Nevertheless, the analysis of the individual spaxels in 
our observations reveals unexpected behaviours when we look 
at the [Nii]/HQf vs. [Oiii]/H/5 diagram, shown in Fig. [S] In order 
to be compared with our observations, we have taken as a refer- 
ence the theoretical curve given by Kewley et al. (2001, in red 
solid line) and the empirical fitting described by Kauff'mann et 



al. (2003), in dashed red line). Following Kewley et al. (2006), 
the region between these two lines could correspond to galax- 
ies with a composite source of ionization, partially from star- 
formation and partially from X-rays, while the area under the 
curves could correspond to pure star-forming regions and the 
area over the curves to pure AGNs. LINERs would lie in the low 
region at right. As can be seen, for all the three objects, most of 
the spaxels in Region 1 (dark blue squares) and the integrated 
spectrum (represented by a large stripped square) lie in the SF 
region, although the points are close to the Kewley 's curve for 
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HS 0837+4717. Nevertheless, a large fraction of the spaxels in 
Regions 2 and 3 (i.e. those in the outer regions of the galaxies) 
lie both in the composite and AGN region in HS 0128+2832 and 
Mrk 930. In the case of HS 0837+4717, most of these spaxels 
lie in the AGN region. 

The excitation of the ionized gas is studied by using the ratio 
of emission-line fluxes of [On] at 3727 A and [Oiii] at 4959, 5007 
A. Although this ratio has a certain dependence on the efl'ective 
equivalent temperature of the ionizing cluster and on the metal- 
licity of the gas, it depends mainly on the ionization parameter 
(i.e. the ratio between the number of ionizing photons and par- 
ticles, Diaz, 1998). In Fig. [6] we see the maps of the logarithm 
of this ratio and their corresponding distributions. The typical 
errors of this ratio range from 0.15 dex in HS 0128+2832 and 
HS 0837+4717 to only 0.03 dex in Mrk 930. This ratio tends 
to be higher as a larger number of spaxels in the outer regions 
is considered, consistently with a larger distance to the ionizing 
central stellar cluster. In the case of Mrk 930 the gaussian fits are 
found in the inner regions. On the contrary, in HS 0837+4717, 
the gaussian fits only the outer regions, although this is possi- 
bly due to a lower number of spaxels (i.e. for all the spaxels in 
the field of view with the required S/N). For the three objects, 
the values in the integrated spectra are lower than the averages 
found in the histograms of the diff'erent regions, especially in HS 
0837+4717. In the case of Mrk 930, almost no diff'erence has 
been found between the excitation of the southern and north- 
ern knots which, besides, are in good agreement with the values 
found for the gaussian fit of the inner-most region. 



Mrk 930 

35 1 — > > > — 




500 1000 1500 2000 

n ([SI I]) (cm"3) 




Fig. 8. Distribution histograms (top) and map (bottom) of the 
electron density from [Sii] in Mrk 930 The solid lines show the 
same contours as in Fig 3. Units are in cm"^. The distribution 
regions are the same defined in Fig.|4] 



3.5. Electron temperature and density 

Due to the low metallicity content of these galaxies, the cool- 
ing of the gas is not very efficient and the electron temperature 
is high. Therefore, the electron temperature of [Oiii] has been 
estimated in a large number of spaxels for our three galaxies 
by means of the flux ratio between the strong nebular emission- 
fines [Oiii] 4959, 5007 A and the auroral weaker fine at 4363 A. 
In order to calculate the electron temperature t([Oiii]), we use 
the expressions from Hagele et al. (2006, 2008) that are based 
on the TEMDEN task under the IRAF nebular package. The dis- 
tribution functions and the maps of this temperature in the three 
observed galaxies are shown in Fig. [7] The average error in the 
individual spaxels are much higher for HS 0128+2832, where it 
reaches 2000 K. For HS 0837+47 17 and Mrk 930 the uncertainty 
is ~ 700 K. A gaussian fit is found for the temperature values of 
all the regions in HS 0837+4717 and in Regions 1 and 2 in Mrk 
930, while no gaussians have been found for any region in HS 
0128+2818, although for this galaxy the mean value of t([Oiii]) 
(13200 K) agrees within the errors with the value measured in 
the integrated spectrum. For HS 0837+4717, an agreement is 
also found between the electron temperature of the integrated 
spectrum (16900 K) and the mean values derived in the three 
studied regions. For Mrk 930 the electron temperatures derived 
in the integrated spectra of both southern and northern knots are 
sfightly lower (12000 and 1 1000 K, respectively), but still within 
the dispersion of the corresponding gaussian. 

We derived the electron densities from the ratio between the 
[Sii] emission lines at 6717,6731 A using the same TEMDEN 
routine. Only for Mrk 930 it has been possible to measure these 
lines with enough confidence in a large fraction of spaxels. In 
the case ofHS 0128+2832, only an upper limit to the density 
can be provided, although a reliable estimate of the electron den- 
sity has been made in the analysis of the integrated spectrum of 
the brightest region of this galaxy, but the errors are too high to 
achieve any precise determination. In the case of HS 0837+4717, 
the redshift of this galaxy placed the [SII] lines outside of the 
spectral range covered in this work. The distribution function 
and the map of the electron density in Mrk 930 are shown in 
Fig.[8l The average error in the spaxels is 50 cm"^. A gaussian 
fit has not been found with enough confidence level for any of 
the analyzed regions. In fact, over the outer southern part of the 
galaxy, the electron density map displays high values, reaching 
up to 1700 cm"^ that is still below the critical density limit for 
collisional de-excitation. This makes the average density of all 
distributions to be higher in the outer regions and the density 
measured in the integrated spectrum of the southern knot to be 
slightly higher than in the northern knots. 

3.6. Oxygen and nitrogen chennical abundances 

In those spaxels with an estimate of the electron temperature it 
is possible to derive chemical abundances following the direct 
method (i.e. the derivation of abundances using the relative emis- 
sion flux of collisional-to-B aimer emission lines and the electron 
temperature). O abundance has been derived by calculating the 
abundances of the main states of ionization in the nebula: O^^ 
and O^, the most abundant under the excitation conditions of 
the gas in these objects. In both cases, we have used expres- 
sions described in Hagele et al. (2006, 2008) based on the task 
IONIC under the IRAF nebular package. O^^ abundances have 
been derived using the emission line intensities of 4959 and 5007 
A relative to llj3 and the electron temperature of [Oiii], which are 
available for a large number of spaxels. 
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Fig. 7. Distribution histograms (top) and maps (bottom) of the electron temperature of [Oiii] in HS 0128+2832, HS 0837+4717, and 
Mrk 930 from left to right respectively. The solid lines show the same contours as in Fig 3. Units are in K. The distribution regions 
are the same defined in Fig.lH 



In the case of abundances, we have used the electron tem- 
perature of [On] and the intensity of [On] 3727 A relative to HyS. 
Since t([Oii]) has not been measured in any of the spaxels, we 
have derived it from the available t([Oiii]) using the relations 
based on photoionization models proposed by Perez-Montero & 
Diaz (2003), which depend on the value of the electron density. 
In those spaxels without an estimate of the electron density (i.e. 
all in HS 0128+2832 and HS 0837+4717 and some of them in 
Region 3 of Mrk 930), we have assumed a value of the electron 
density of 100 cm"^. The distribution functions and the maps 
of the total O abundance as derived using the direct method are 
shown in Fig. [9l As the main source of uncertainty in this to- 
tal abundance comes from the temperature, the average error is 
higher in HS 0128+2832 (0.3 dex) than in HS 0837+4717 and 
Mrk 930 (0.09 and 0.08 dex respectively). As in the case of elec- 
tron temperature a gaussian fit is found in the three regions of 
HS 0837+4717 and in Regions 1 and 2 of Mrk 930, but in none 
of HS 0128+2832. The mean value of these distributions agree 
within the Icr dispersion range with the values derived from the 
integrated spectrum, although in the case of Mrk 930 there is a 
slight off'set between the O abundance of the southern and north- 
ern knots, but it is within the observational errors. 

A deeper study of the metallicity distribution of these galax- 
ies has been performed using strong-line methods. These do not 
depend on the previous determination of the electron tempera- 
ture and they can be applied in a larger number of spaxels. In the 
case of O abundance there are several indicators based on the 
intensity of the brightest lines of the spectrum. In our case, in or- 
der to be consistent with our previous estimation of the O abun- 
dance based on the direct method, we have chosen the calibra- 
tion proposed by Pilyugin & Thuan (2005), which is calibrated 
using objects whose metallicities have been derived using elec- 



tron temperature estimates. This method relies on the calibration 
of the parameters R23 and P, depending on [On], [Oiii], and H/? 
emission-line fluxes. One of the main drawbacks of the strong- 
line methods based on R23 is the degeneracy of the parameter 
with O abundance (i.e. a same value of R23 leads to two difl'erent 
values of the O abundance). Nevertheless, since our three galax- 
ies are metal-poor objects, the calibration proposed by Pilyugin 
& Thuan for the low metallicity regime has been used. In Fig. 
[TOl we show the distribution functions and the maps of the total 
O abundance in the three studied galaxies using this method. In 
all cases, the average errors in the spaxels of the three galaxies 
are lower than the uncertainty associated with the calibration of 
this parameter, which is 0.2 dex according to Perez-Montero & 
Diaz (2005). As in the case of the O abundance derived using 
the direct method, in all three galaxies difl'erent gaussian func- 
tions are fitted at diff'erent scales. In the case of the inner-most 
region of HS 0837+4717, the level of significance is very large 
(98%), but it decreases substantially for the outer regions. In the 
case of HS 0128+2832, a gaussian fit is found only in Region 2. 
Finally, in Mrk 930, contrary to the direct method where gaus- 
sians were found in Regions 1 and 2, a normal function is only 
fitted in Region 1. For the three objects, a good agreement is 
found between the mean values of these distributions and their 
equivalents using the direct method in a lower number of spaxels 
and with the metallicity estimates from the integrated spectrum 
using the same parameters. 

Nitrogen abundances have been also calculated using both 
the direct method in those spaxels with an estimate of the elec- 
tron temperature and using a strong-line method for all the other 
spaxels where a good measure of the bright emission lines ex- 
ist. In the first case, we have used the relative intensity of the 
[Nil] 6584 A to UJ3 and the electron temperature of [Nii], esti- 



10 



E. Perez-Montero et al: A PMAS study of high N/0 BCDs 




Mrk 930- 12+log(O/H) 





I 




Fig. 9. Distribution histograms (top) and maps (bottom) of the total oxygen abundance, in terms of 12+log(0/H) calculated using 
the direct method, in HS 0128+2832, HS 0837+4717, and Mrk 930 from left to right, respectively. The solid lines show the same 
contours as in Fig. 3. Units are in dex. The distribution regions are the same defined in Fig.lH 



mated using the relation with t([Oiii]) based on models described 
in Perez-Montero & Contini (2009). The N/0 abundance ratio 
is then assumed to be equal to the N^/0^ ratio. The distribu- 
tion functions and the maps of this ratio are shown in Fig. [TT] 
As in the case of O total abundances derived from the direct 
method, the mean error of the N/0 ratio in each spaxel is higher 
in HS 0128+2832 (0.4 dex), than in HS 0837+4717 (0.2 dex) 
and Mrk 930 (0.16 dex). Contrary to all the other distributions 
in this work, gaussian functions fit with a good level of confi- 
dence in all the regions of the three galaxies, in agreement with 
the values derived in the integrated spectra. 



For the rest of spaxels, with no derivation of the temperature, 
the N/0 ratio can be estimated by resorting to the empirical cal- 
ibration of the N202 parameter, proposed by Perez-Montero & 
Contini (2009). This parameter is based on the ratio of [Nii] and 
[On] emission lines, presenting a linear relation with the N/0 
ratio. As previously, we have used this method because it has 
been calibrated using data with a direct estimation of the N/0 
ratio and it is consequently consistent with the N/0 derivation 
from the electron temperature. The distribution functions and the 
maps of the N/0 values using the strong-line method are shown 
in Fig. [121 The average errors are in all the spaxels lower than 
the uncertainty associated with the calibration of this parameter, 
which is 0.2 dex, according to Perez-Montero & Contini (2009). 
The mean values of the N/0 ratio distributions as derived from 
this strong-line method are compatible with the values found 
using the direct method. Nevertheless, the gaussian fits appear 
only in the inner regions of HS 0837+4717 and no gaussians are 
found in any region for Mrk 930 and HS 0128+2818. 




Fig. 13. Radial velocity of the ionized gas for the Ha emission 
line in Mrk 930. The solid lines show the same contours as in 
Fig. 3. Units are in km/s. 



3. 7. Kinematics of Mrk 930 

We have studied the velocity field of the ionized gas by fitting a 
gaussian to the Ha emission line in the three galaxies, but only 
Mrk 930 presents a rotation curve, as shown in the radial velocity 
map of Fig. \T3\ The spectral resolution of our data prevents us 
from deriving reliable velocity dispersions. The west side of the 
galaxy appears to be less disturbed than the east side; a sharp 
velocity gradient in the north- south direction is visible in the 
eastern half of the map. This velocity field could be the result of 
an interaction system. 
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Fig. 10. Distribution histograms (top) and maps (bottom) of 12+log(0/H) calculated using the strong-line method based on the 
Pilyugin & Thuan calibration, in HS 0128+2832, HS 0837+4717, and Mrk 930 from left to right, respectively. The solid lines show 
the same contours as in Fig. 3. Units are in dex. The distribution regions are the same defined in Fig.|4] 



4. Discussion 

4.1. On the spatial uniformity of tlie derived properties in tlie 
sample 

We have studied the distribution of diff'erent physical properties 
and chemical abundances along the field of view of the PMAS 
instrument in some BCDs with diff'erent degrees of excess in 
their N/0 ratio in relation to the expected values for their metal- 
licity regime. IPS allows a diff'erent approach in their study in 
contrast to integrated long-slit or fiber observations. In these 
cases, the integrated spectrum is a weighted mean of the dif- 
ferent regions of the observed objects, whose statistical weights 
are the relative luminosities of the diff'erent regions. On the con- 
trary, IPS data can be analyzed considering each observed region 
independent of their relative brightness. In this way, in the dis- 
tributions of the physical properties and chemical abundances of 
the three studied BCDs in this work, it is remarkable the diff'er- 
ences between the mean of the distributions in diff'erent regions 
and the values of the same quantities measured for the integrated 
spectrum. These diff'erences are well illustrated by the redden- 
ing constant distribution in the HS 0128+2832 galaxy (see first 
panel of Fig.|4l), whose integrated value is c(HyS) = 0.14, while 
the mean value of the spaxels distribution is c(Uj3) = 0.37 in the 
inner-most region. This diff'erence is understood by inspecting 
the shape of the reddening constant distribution, where no gaus- 
sian fit is found for any of the regions and, on the contrary, an 
apparent double-peak distribution appears, with one of the max- 
ima peaking at the same value found for the integrated spectrum 
and the other maximum peaking at the same value that the mean 
of all the data. Another example of this disagreement is found in 
the study of the density distribution of Mrk 930, which is far to 
be fitted by a gaussian. This is partially motivated by the pres- 



ence of a high density structure well visible in the south part of 
the galaxy and which is also seen in the histogram distribution 
plot (upper panel of Fig. (8]). This structure is also responsible for 
the diff'erence found in the electron density between the north 
(32 cm"^) and south (< 190 cm"^) parts of the galaxy. Besides, 
the mean values found for diff'erent regions of the same galaxy 
range between 101 and 271 cm"^, as we go from Region 1 (i.e. 
the inner-most region, around the peak in Ha) and Region 3 (i.e. 
all valid spaxels in the IFU), as we can see in Tabled In the most 
part of the cases, this disagreement between the value found in 
the integrated spectrum and the mean value of the spaxels distri- 
bution appears when no gaussian fit is found for this distribution, 
indicating local spatial variations which can not be understood 
only in terms of statistical fluctuations. 

The histograms of the properties studied in difl'erent regions 
have been also used to find out to what extent a certain prop- 
erty can be considered as homogeneous within the observational 
and statistical errors in diff'erent regions. This would allow us to 
find in principle possible variations of the excitation as a func- 
tion of the distance to the ionizing cluster, to find regions with 
excesses in electron density or extinction, or to locate the chem- 
ical pollution in specified places of the observed field. However, 
the assumption of the hypothesis that a certain parameter can 
be considered as uniform if its spaxels distribution can be fitted 
by a gaussian have two important limitations in this work which 
must be taken into account. Firstly, there is a sampling eff'ect 
which makes that, for some very compact objects, the number 
of spaxels to be considered as valid can be somehow reduced. 
This is the case of the distributions of excitation and N/0 ratio, 
as derived using the N202 parameter, in HS 0837+4717 (middle 
panels of Fig. [5] and [121 respectively), in which the null hypoth- 
esis is rejected for Region 1. However, a gaussian fit is found 
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Fig. 11. Distribution histograms (top) and maps (bottom) of log(N/0) calculated using the direct method, in HS 0128+2832, HS 
0837+4717, and Mrk 930 from left to right, respectively. The solid lines show the same contours as in Fig. 3. Units are in dex. The 
distribution regions are the same defined in Fig. |4] 



when the number of spaxels grow and Regions 2 and 3 are also 
considered. This is possibly indicating that the rejection of the 
gaussian fit hypothesis does not imply that a certain parameter is 
not uniform in a certain region if there is a low number of points 
to be considered in the corresponding statistics. The second lim- 
itation comes from the poor seeing measured during the night of 
observations, which is worse in all cases than the working spa- 
tial resolution. In this case, serious doubts appear about the no 
correlation of the fluxes between adjacent spaxels. Nevertheless, 
the typical size of the regions where the uniformity has been as- 
sumed is much larger in all cases than the typical value of the 
point spread function for the quoted seeing, and this is indicative 
of the validity of the results. 

Another important fact which reinforces the assumption of 
uniformity in those cases where a gaussian can be fitted, is that 
most of the times the normal functions are found in more than 
a region in each galaxy. With the exception of the O abundance 
as derived from the R23-P method in HS 0128+2832 and the 
N/0 ratio as derived from the N202 method in Mrk 930, the rest 
of the physical properties and chemical abundances are fitted in 
two or even in the three defined regions in the field of view. This 
excludes the possibility in those cases of an artificial gaussian 
distribution in a defined region of the space, as the fitting is found 
simultaneously in diff'erent and independent regions. 

The results of these fittings listed in Table 3 indicate that in 
HS 0128+2832 only the N/0 ratio as derived following the direct 
method, can be fitted by a normal distribution. On the contrary, 
in HS 0837+4717, all the distributions can be fitted by a normal 
function at diff'erent regions, with the exception of the excitation 
and N/0 ratio as derived from the N202 method in Region 1, as 
explained above. Finally, in Mrk 930 only reddening and density 
distributions are not fitted by a gaussian. In this case, it is the 



disturbed morphology of this galaxy the most probable cause of 
the found variations. 



It is remarkable the fact that in all the three galaxies the N/0 
ratio as derived from the direct method can be considered as uni- 
form at the spatial scales where it can be measured with enough 
precision, independently of the mean value of their distribution 
or the possible presence of WR stars. However, this is not the 
case for the same ratio when it is derived using an empirical cal- 
ibrator, N202 in this case, which does not appear as uniform 
in several regions of HS 0128+2832 and Mrk 930. For total O 
abundances, the same discrepancy between the distributions de- 
rived from the direct method and from the empirical calibration 
of the R23 parameter given by Pilyugin & Thuan (2005) is ob- 
served. Hence, in the case of Mrk 930, the O abundance distribu- 
tion is well fitted by a gaussian when it is derived using the direct 
method in Regions 1 and 2, but the gaussian is only fitted in the 
inner-most region when this abundance is derived using R23. 
Several limitations related to the calibration of these parameters 
can be taken into account to explain this discrepancy, including 
the sample used to calibrate them, mostly data coming from inte- 
grated observations of star forming knots, and the dependence of 
the calibrators on other functional parameters, such as ionization 
parameter or eff'ective equivalent temperature (Perez-Montero & 
Diaz, 2005), which are not uniform along a galaxy. Therefore, 
although the strong-line parameters allow the sampling of chem- 
ical abundances in a larger number of positions, they are less 
suitable to study the homogeneity of these abundances in com- 
parison to other direct method-based techniques. 
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HS 0837+4717 




Fig. 12. Distribution histograms (top) and maps (bottom) of log(N/0) calculated using the strong-line method based on the ratio 
[NII]/[OII], in HS 0128+2832, HS 0837+4717, and Mrk 930 from left to right, respectively. The solid lines show the same contours 
as in Fig. 3. Units are in dex. The distribution regions are the same defined in Fig.lH 



4.2. Excitation variations and their Implications on diagnostic 
diagrams 

One of the most interesting results extracted from the analysis 
of the spatial distributions of the properties of these three stud- 
ied objects appears with the position of several spaxels in Fig. [51 
which represents the [Nii]/Hck vs. [Oiii]/HyS diagnostic diagram. 
As we see, the position of the integrated points and the spaxels 
belonging to Region 1 lie in the star forming region, although the 
proximity of these points to both the theoretical and empirical 
curves are sensitive to the N/0 ratio (i.e. in HS 0128+2832 and 
HS 0837+4717, which present the largest relative N abundance, 
these points are closer than in Mrk 930). Nevertheless, a frac- 
tion of the outer spaxels, which is larger in HS 0837+4717 and 
lower in Mrk 930, spread towards right in the diagrams so they 
can lie in the so-called composite region or even in the active 
galactic nuclei region. This eff'ect is probably due to the com- 
bination of two diff'erent eff'ects: i) the high relative N/0 ratio 
which makes some points to lie in these regions, as shown al- 
ready by Perez-Montero & Contini (2009), and ii) the excitation 
structure of the galaxy, in which the [Nii] emission line fluxes 
remain bright even at distances where the other emission lines 
are sensibly less bright. This makes the [NiiJ/Hof to be higher in 
these regions, while [Oiii]/H/5 remains constant or slightly de- 
creases. In this case the spectra of some positions of the field of 
view could be wrongly classified as AGNs. This eff'ect was al- 
ready observed by James et al. (2010) in the BCD Mrk 996 using 
also IFU observations. 



4.3. N/0 excess and the presence of WR stars 

Although the WR blue bump has been reported in the integrated 
spectra of HS 0837+4717 (Pustilnik et al., 2004) and Mrk 930 



(Guseva et al., 2000), we have not detected it neither in any of the 
individual spaxels of our IFU observations nor in the integrated 
spectra of the brightest regions. The low S/N of the bump com- 
bined with poor seeing conditions during the observing night 
could be the main reasons of this fact. Nonetheless, it is possible 
to carry on the spatial analysis of the physical properties and the 
O and N chemical abundances along the field of view of the IFU. 
The analysis described in section 3.6 shows that, in most of the 
regions of the three galaxies, we can fit a gaussian function to the 
O abundance and N/0 ratio distributions, obtained from both the 
direct method and the strong-line calibration from R23-P for O 
abundance and N202 for N/0 ratio. Besides, the mean values of 
these gaussian functions match within the errors the values de- 
rived using the integrated spectra from both the brightest region 
and the reported long-slit observations in each galaxy. These val- 
ues are consistent with a scenario of low metallicity with an ex- 
cess of the N/0 ratio in relation to the expected value at this 
metallicity regime. Adopting log(N/0) = -1.6 as the most typ- 
ical value of the plateau in the 12+log/O/H) vs. log(N/0) rela- 
tion, as can be seen in Fig. 1, this excess is Alog(N/0) = +0.43 
dex in HS 0128+2832, +0.78 dex in HS 0837+4717 and +0.17 
dex in Mrk 930, taking as a reference the mean value of the 
gaussians fitted to the distributions of log(N/0) as derived us- 
ing the direct method in the respective Region 3 of the galaxies. 
At same time, we can take into account the mean radii of the 
spatial distributions of log(N/0), considering the spatial size of 
the spaxels at the adopted distances for each galaxy, which are 
330 pc in HS 0128+2832, 840 pc in HS 0837+4717, and 380 
pc in Mrk 930. By inspecting the corresponding maps, the sizes 
of the log(N/0) distributions as derived from the direct method 
are approximately 1.2 x 1.7 kpc^ for HS 0128+2832, 2.5 x 3.4 
kpc^ for HS 0837+4717, and 1.9 x 3.0 kpc^ for Mrk 930, assum- 
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Fig. 14. Evolution of the N/0 ratio for a spherical gaseous distribution of constant density with an initial metallicity of 12+log(0/H) 
= 8.0 and log(N/0) = -1.6, which is polluted by the stellar winds coming from a stellar cluster with Z = 0.004. At left, we show the 
effects of a cluster of mass log(M/Mo) = 7 for different radii of the distribution. At right, we show the effects of clusters of different 
stellar masses on a gaseous distribution with a radius of 1 kpc. Stellar wind yields have been taken from Molla & Terlevich (201 1, 
in prep.) 



ing that these are upper limits due to the seeing of the observing 
night. 

These scale lengths lead us to question if the reported WR 
stars can be the actual source responsible for the global N pol- 
lution detected in our sample of BCDs. To check the validity of 
this possibility we have taken the chemical yields of different 
elemental species ejected to the ISM by massive stellar clusters 
proposed by Molla & Terlevich (in prep.) for a metallicity of the 
cluster Z = 0.004 (^ l/SZ©), which is close to the metallicity of 
our 3 observed BCDs. We have studied the effects of these winds 
on a spherical homogeneous distribution of gas with constant 
density of 100 cm"^ and an initial O abundance of 12-Flog(0/H) 
= 8.0 and log(N/0) = -1.6. We also assume that the mixing 
of the ejected elements with the material of the spherical distri- 
bution is instantaneous, which is not reallistic at these galactic 
scales, but it is the most favourable hypothesis for the assump- 
tion of the pollution of the ISM by stellar winds. 

The evolution of the log(N/0) ratio with time is shown in 
Fig. [14] under different assumptions. In the left panel, we show 
the effect of the ejecta from a massive cluster with log(M/Mo) = 
7 to a spherical gaseous distribution of different radii. As we can 
see, for the typical radii of the N/0 distributions in our galaxies, 
which are of the order of 1 kpc, the effect is negligible. For lower 
radii, it begins to be apparent how the ejection of this material 
changes the relative chemical composition of the ISM. At first, 
the ejection of O by very massive stars makes the N/0 ratio to 
decrease and the later appearance of WR star winds makes the 
N/0 ratio to increase again, reaching a maximum and an asymp- 
totic value. The value of this maximum matches only the mea- 
sured mean values of our three galaxies for radii much lower 
than those measured in the IFU observations. In the case of HS 
0837-1-47 17, the measured N/0 is even higher than the maximum 
value of log(N/0) reached at any radius. 

In the right panel of Fig. [14] we have fixed the radius of the 
spherical gas distribution to 1 kpc, and we have varied the stellar 
mass of the ionizing cluster. As can be seen, the effects are only 
noticeable at this scale when the cluster is more massive than 
log(M/Mo) = 8. We can estimate the masses of the ionizing stel- 



lar clusters in the three BCDs, taking the following expression 
proposed by Diaz (1998) 

log(M/Mo) = logL(Ha) - 0.86 • EW(Hp) - 32.61 

with the Ha luminosities measured from the area covered by the 
IFU and the EW(I^) of the integrated spectrum of the corre- 
sponding brightest regions. Therefore, the ionizing masses are 
of 106-47 Mo for HS 0128-h2832, lO^ ^^ M© for HS 0837-h4717, 
and lO^-^^ M© for Mrk 930. Although these estimates are only 
lower limits, because escaping photons and dust absorption are 
not taken into account, they are much lower than the mass of the 
ionizing cluster required to pollute the ISM to the spatial scales 
measured in our IFU observations. 

Therefore, although a number of WR stars have been de- 
tected in two of these galaxies (i.e. of the order of 10^ both in HS 
0837-h4717 by Pustilnik et al, 2004 and in Mrk 930 by Guseva 
et al., 2000), and according to stellar evolutionary models they 
are also expected to be present in HS 0128 -h 28 3 2, the global high 
N/0 ratio observed in these objects is not likely produced by the 
pollution from stellar wind ejecta coming from WR stars, even 
in the most favourable conditions of mixing and dillution and, 
on the contrary, could be more related to other global processes 
affecting the metal content of the whole galaxy, as was already 
suggested for the green pea galaxies analyzed by Amorm et al. 
(2010). They propose a combination of outflows of enriched gas 
and inflows of metal-poor gas, capable to trigger the star for- 
mation processes in these objects, decreasing the total content 
of metals without noticeably changing the N/0 ratio (Koppen & 
Hensler, 2005). 



5. Conclusions 

We have studied, by means of PMAS optical integral field 
spectroscopy, the spatial distribution of ISM physical-chemical 
properties for a sample of three BCDs (HS 0128-h2832, HS 
0837-1-4717, Mrk 930) with low metallicity and an excess of 
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N/0, as measured with previous integrated long-slit spectropho- 
tometry. The spectral range covering from 3700 to 6900 A and 
the spatial resolution of l"xl" in a field of view of 16"xl6" en- 
sure the characterization of the physical properties and the chem- 
ical abundances across the objects, although the spatial resolu- 
tion is somehow reduced due to the seeing during the observing 
night 2"). 

Ua intensity maps show that both HS 0128-h2832 and HS 
0837-^4717 have compact aspect, with their luminosities domi- 
nated by a single central massive stellar cluster. Mrk 930, whose 
morphology has been also compared with available UV HST- 
ACS images, shows a more complex structure, with knots in 
the northern and southern part of the galaxy. Besides, it shows a 
clear rotation curve along the north-south axis, which could be 
consequence of an interaction system. 

All properties derived from the emission-line intensities with 
enough S/N have been analyzed in diff'erent spaxels of the field 
of view. Their corresponding distribution in diff'erent regions 
have been also compared to the properties of the objects as de- 
rived from the integrated spectrum collected in the brightest re- 
gion of each galaxy. For instance, while all integrated spectra lie 
in the star-forming region of the [Nii]/Ha vs. [Oiii]/HyS diagram, 
a sample of the spectra in the outer regions of each galaxy ex- 
tends towards right in this diagram, even lying in the so-called 
composite and AGN regions. This is due to the high N/0 ratio of 
these spaxels, combined with the excitation structure of the gas, 
causing the [Nii] emission lines to remain bright even at very 
large distances to the ionizing stellar source. 

We have assumed that a certain property of the ISM can be 
considered as uniform within the galaxy or a determined region 
if the null hypothese of the Lilliefors test (i.e. the fitting of a 
gaussian function to the corresponding distribution) cannot be 
rejected. We find that in those cases where a gaussian does not fit 
a given distribution, the disagreement between the mean value of 
the distribution and the value derived from the integrated spec- 
trum is larger. 

We find that extinction, as derived from the decrement of 
Balmer, is only uniform in HS 0837-h4717. Mrk 930 shows a 
substructure of higher extinction in the outer southern part of 
the galaxy which breaks the uniformity. Excitation, as estimated 
using the [Oii]/[Oiii] emission-line ratio is only uniform in the 
inner-most regions of HS 0837-h4717 and Mrk 930, due possi- 
bly to the lower excitation in the outer parts of these galaxies, 
as the distance to the corresponding ionizing sources increases. 
Electron density has been only estimated with enough confi- 
dence in Mrk 930, where a structure of high density appears in 
the south part of the galaxy, so the distribution of this property 
cannot be considered as uniform. 

Electron temperature has been derived in a sample of spaxels 
of the three galaxies from the [Oiii] 4363 A to [Oiii] 4959,5007 
A emission line-ratio. Its distribution appears as uniform in HS 
0837-h4717 and Mrk 930. This is also the case for O chem- 
ical abundances and N/0 ratios as derived using these elec- 
tron temperatures through the direct method. In the case of HS 
0128-1-2832, the N/0 derived using the direct method is the only 
ISM property in this galaxy which is uniform after our analy- 
sis. The distribution of 0/H and N/0 derived using strong-line 
methods are uniform in a lower number of regions, possibly due 
to the dependence of the corresponding strong-line parameters 
on the ionization structure of the gas. 

Finally, we have not detected neither the broad WR blue 
bump (~ 4650 A) nor the red one (~ 5808 A) suggested by previ- 
ous works in two of the observed galaxies. Even though we have 



investigated if these stars can be considered as the main cause of 
the measured excess of N/0. We have demonstrated, using the 
stellar yields of massive star winds at the corresponding metal- 
licity of these objects (Z=0.004, Molla & Terlevich, in prep.) 
and taking into account the stellar masses of the ionizing clus- 
ters as estimated from the integrated Ha flux, that at the scale 
lengths where we detect the uniform N/0 high values (i.e. ~ 1 
kpc), the WR stars cannot be responsible for the enhancement in 
the N abundance. It would be necesary between 2 and 3 orders 
of magnitude higher in the ionizing cluster mass to reach such 
a degree of pollution, assuming an instantaneous mixing of the 
ejected material with the surrounding gas. Thus, for the case of 
these three studied galaxies, another chemical evolution scenario 
(e.g. metal-rich outflows, infall of metal-poor gas or interacting 
scenario as the morphology and dynamics of Mrk 930 suggest) 
is required. 
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